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The potential and current distributions are 
calculated across the width of a quantum 
Hall effect sample for applied currents be- 
tween (xA and 225 (jlA. For the first 
time, both a confining potential and a cur- 
rent-induced charge-redistribution potential 
are used. The confining potential has a 
parabolic shape, and the charge-redistribu- 
tion potential is logarithmic. The solution 
for the sum of the two types of potentials 
is unique at each current, with no free 
parameters. For example, the charge-deple- 
tion width of the confining potential is de- 
termined from a localization experiment by 
Choi, Tsui, and Alavi, and the spatial extent 
of the conducting two-dimensional electron 
gas across the sample width is obtained 
from the maximum electric field deduced 
from a high-current breakdown experiment 
by Cage and Lavine, and from the quantum 



Hall voltage. The spatial extent has realistic 
cut-off values at the sample sides; e.g., no 
current flows within 55 magnetic lengths of 
the sides for currents less than 215 \lA. 
The calculated potential distributions are in 
excellent agreement with contactless elec- 
tro-optic effect laser beam measurements of 
Fontein et al. 



Key words: breakdown of dissipationless 
state; charge distributions; current distribu- 
tions; logarithmic charge-redistribution po- 
tential; parabolic confining potential; poten- 
tial distributions; quantum Hall effect; 
two-dimensional electron gas. 



Accepted: June 20, 1995 



1. Introduction 

The potential and current distributions within quantum 
Hall samples are important aspects of the integer quan- 
tum Hall effect [1-3]. In this effect the Hall resistance, 
Rh, of the ith plateau of a fully quantized two-dimen- 
sional electron gas (2DEG) has the value Riiii)=hlie^i), 
where h is the Planck constant, e is the elementary 
charge, and i is an integer. Early attempts to measure 
potential distributions across samples [4-7] used electri- 
cal contacts to the two-dimensional gas that were placed 
within the sample interior. The potentials were found to 
vary throughout the entire sample. There was concern, 
however, that the electrical contacts themselves signifi- 
cantly altered the potential distributions. Fontein et al. 
[8] have made contactless measurements of potential 



distributions using a laser beam and the electro-optic 
Pockels effect. They observed major fractions of the 
quantum Hall voltage occurring near the sides of the 
sample, but also significant contributions within the 
interior. Valid predictions of the potential distribution 
across quantum Hall samples should agree with their 
results. 
In this paper we calculate the potential distributions 
across the sample for applied currents /jd between (xA 
and 225 (xA by: (a) assuming a parabolic confining 
potential for the charge carriers and using parameters of 
the parabola obtained experimentally by Choi, Tsui, and 
Alavi [9] ; (b) assuming an applied current-induced loga- 
rithmic charge-redistribution potential for the charge 
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carriers that is similar to that of Balaban, Meirav, and 
Shtrikman [10], but with very different cut-off values 
for the spatial extent of the potential; (c) assuming that 
the width of the conducting region varies with applied 
current because a Lorentz force deflects the conducting 
electrons slightly towards one side of the sample; (d) 
using the maximum electric field deduced by Cage and 
Lavine [11] from a breakdown experiment at high cur- 
rents to determine the cut-off value for the spatial extent 
on one side of the sample; and (e) using the quantum 
Hall voltage to determine the cut-off value of the spatial 
extent on the other side of the sample. The calculated 
potential distributions are in excellent agreement with 
the measurements of Fontein et al. [8]. 



2. Coordinate System 

The coordinate system is shown in Fig. 1 . For conve- 
nience in writing the equations, the origin is located at 
the source S and is halfway across the sample width w. 
The sample labeling is chosen to be consistent with 
previous work [11]. Potential probes 1 and 2, and the 
drain D, are not shown. The positive x axis points along 
the sample in the general direction of the externally 
applied current /so. The positive y axis is chosen as 
indicated. Therefore the positive z axis points downward 
for a right-handed coordinate system. The magnetic 
field B also points downward, simply to be consistent 
with results from the breakdown experiment [11] that 
will be used in Sees. 4 and 5. 

Note that the conducting charges are electrons with 
charge q = - e. This is taken into account throughout the 
paper; it is necessary to do so because the signs of both 
the confining potential and the charge-redistribution 
potential depend on the sign of the charge carriers. 



W/2 




F^A > — 



E/^Z/ 



Fig. 1. The coordinate system used in this paper. Only part of the 
sample is shown. The origin is located at the source S, halfway across 
the sample width w. The dotted curves indicate the electron flow 
pattern for this magnetic field direction, f l is the Lorentz force on the 
conducting electrons and Fc is the Coulomb force. B is the magnetic 
field, V is the electron velocity, and Va and Vb are the potentials on 
either side of the sample. 



The Lorentz force Fl = ev X B is in the positive y 
direction. This force deflects the electrons slightly to the 
right until it is matched by the Coulomb repulsive force 
Fc = -eE [12]. A charge-redistribution of the 2DEG 
results from this deflection. Also because of the Lorentz 
and Coulomb forces, the electrons enter the corner of the 
source aty = -w/2 for this magnetic field direction and 
exit at the corner +w/2 of the drain — in agreement with 
the experiment of Klass et al. [13]. We assume that the 
electrons spread out across the sample interior in agree- 
ment with the experiment of Fontein et al. [8]. Potential 
probes 4 and 6 are near the potential of the source. 
Probes 3 and 5 are near the potential of the drain, and 
have a positive potential relative to the source for these 
current and magnetic field directions. The chemical po- 
tential (^A = Va is therefore positive relative to the chem- 
ical potential 0b = Vb on the opposite side of the sample. 



3. Confining Potential 

We begin the calculations with a confining potential 
to prevent the 2DEG from spilling out the sides of sam- 
ples. Choi, Tsui, and Alavi [9] performed an experiment 
on mesa-etched GaAs/AlGaAs hetero structure samples 
in zero magnetic field. They then used one-dimensional 
localization theory to deduce the charge-depletion 
widths, A , of the confining potentials, and found that A 
was (0.5 ± 0.2) [jum for a 2DEG of surface number 
density n^ = i(eB/h) = 1.5 X lO'Vcml We will use their 
results to define the depletion width of the confining 
potential for a mesa-etched sample. 

3.1 Charge-Depletion Region 

Figure 2 (a) shows a schematic of the charge distribu- 
tion in the GaAs/AlGaAs interface region near one side 
of the mesa when there is no applied magnetic field. The 
GaAs layer of our sample [11] has a residual donor 
density of about 1 X lO'^'/cm^, while the donor concen- 
tration in the AlGaAs layer is about 1 X lO'Vcm^ and 
Ms = 5.94 X lO'Vcm^. There is an ionized donor atom in 
the AlGaAs layer for every electron in the 2DEG but, 
unlike Choi, Tsui, and Alavi [9], we assume the ionized 
donor atoms are distributed over a volume rather than in 
a surface sheet with density n^. The confining potential 
is generated from electron surface charges on the side of 
the mesa, as indicated in the figure. There is an ionized 
donor atom or ionized impurity site in the charge-deple- 
tion region for every surface charge. 

We assume a homogeneous charge-depletion region 
in Fig. 2 (b). The depletion width A for a homogeneous 
three-dimensional material is [14] 



A = (2e,VJeNo)' 



(1) 
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Fig. 2. (a) Schematic diagram of the GaAs/AlGaAs interface region 
near one side of the mesa. See Sec. 3. 1 for further explanation, (b) The 
ionized donor charge-depletion density distribution No- (c) The con- 
fining potential V^ for negatively charged surface states. 



where e^ = keo is the dielectric permittivity of the semi- 
conductor, K is the dielectric constant (k = 13.1 for 
GaAs), Eo is the permittivity of vacuum, Vm is the value 
of the confining potential at ±w/2 and A^d is the average 
density of ionized donors and impurity sites in the 
charge-depletion region. We selected the value of the 
charge-depletion width to be Zl = 0.5 (xm [15]. This 
value is consistent with the results of Choi, Tsui, and 
Alavi [9] . We chose the value of Vm to be one-half the 
1.50 V separation between the valence and conduction 
bands of GaAs at 1 K [14], or V^ = 0.75 V. The value of 
the average charge-depletion density from Eq. (1) is thus 
A^D = 4.3 X lO^Vcm^, which seems quite reasonable. 

3.2 Confining Potential Equation 

A homogeneous charge-depletion region results in a 
parabolic confining potential Vc, with the origin at 
J = A = w/2-A, as indicated schematically in Fig. 2 (c). 
The confining potential is negative because the charges 
on the side of the mesa are electrons. 

The equations for the confining potential Vc and its 
electric field E^ = - W. are 



for A < y : 



KCy) = and E,(y) = 



for - A < J < A , 



(2b) 



V,(y) = - a(y+Xf and E,(y) = 2a(y+X) (2c) 



for-— <-y <-A, 



where a = VJA^ = 3.0 X 10" V/m" for A = 0.5 |xm and 
Vm = 0.75 V, and 



A = I - 4 . (3) 



3.3 Confining Potential at /sd = /tt A 

Given the values of A and V^, there is a surprising 
amount that can be deduced about the electron states of 
the confining potential when the magnetic field is ad- 
justed to be at the center of the / = 2 quantum Hall 
plateau and /sd = (xA. Since there is no applied current, 
and therefore no Hall voltage, the Fermi energy Sp is 
constant across the sample width and is located halfway 
between Landau levels. Under these conditions, states of 
the lowest (A^ = 0) Landau level are occupied up to the 
Fermi energy e^ = AcoJ2, no states are occupied in the 
second (A^ = 1 ) Landau level, co^ = eB/m * is the 
cyclotron angular frequency, m * is the reduced mass of 
the electron (0.068 times the free electron mass in 
GaAs), and ^ = h/2TT. References [11,16-18] describe 
how these states can be defined in the Landau gauge. 

Figure 3 shows a schematic drawing of the energy of 
the confining potential for /sd = [xA, with greatly 
exaggerated values of A and ^co^, and only a small frac- 
tion of allowed states. The occupied/unoccupied states 
are indicated as solid/open circles, and the occupied 
(filled) states lie between j^ax and y^m = -Jmax- In the 
presence of the magnetic field, electrons of the 2DEG 
occupy Landau level states that penetrate into the 
charge-depletion regions near the mesa edge, and cur- 
rent circulates around the sample periphery. Under these 
conditions 



_ ^ _ 



eKCymax) = eaCymax-A)^ 



Vn . 



A)^ 



(4) 
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where A = w/2 - A. 

The occupied states of the right-hand side (rhs) con- 
fining potential generate a total current /^ (rhs) that is 



Jmax Jma: 

J J^(y) dy= j 



/c(rhs) = J^(y)dy= (T^E^(y)dy. 



■ i- [KO'n.J - K(A)], = - %=) ^ (5) 



where Jc(y) is the current density, a^y is the off-diagonal 
conductivity tensor component, KCymax) = -aCVmax- 
A)^, and Vc(A) = 0. In the absence of significant dissipa- 
tive scattering on the quantum Hall plateau, cr^y = 1/Rh 
[12]. Similarly, 



-A -A 

/c (Ihs) = Jcdy = cr^E^iy) dy = 



■ j- [K(-A) - KO'n.in)], = - %^ , (6) 



where Vc(ymm) = -a(yMn+y^f- 

It follows from Eqs. (2) to (6) for the 12 906.4 D,, 
i = 2 plateau at 12.3 T, for the 400 (xm wide sample of 
Ref. [11], and for ho = (xA that 



, ,, , /io}, ie^B ^„, , , „, , ,^, 

'^ (*^> = 2^ = 4^ = O-^l f^^ = - ^= ^'""'^ ' ^'^ 




y (urn) 



Fig. 3. Schematic drawing of the energy of the confining potential 
Vc across the sample when ha = (jlA. Values of the charge-depletion 
width A and the Landau energy level spacing Aia^ are greatly exagger- 
ated. The occupied/unoccupied states of the first two Landau levels are 
shown as solid/open circles. The occupied (filled) states lie between 
the locations y^^^ = — j^n. 



4. Charge-Redistribution Potential 

Section 2 noted that the Lorentz force exerted on the 
conducting electrons causes deviations - eha(y) from 
the average surface charge density - ea^^^ = -en^ = 
- ie^Blh of the 2DEG charge-redistribution across the 
sample width. The resulting charge-redistribution 
potential, V^iy), arising from applied currents would be 
a linear function of y if the mobile electrons occupied a 
three-dimensional volume. They occupy a two-dimen- 
sional sheet, however, and MacDonald, Rice, and 
Brinkman [19] expressed this charge-redistribution self- 
consistently in terms of a charge-redistribution potential 
as 



Jmax = - ymm = 199.559 |Jtm , (8) 



and 



Vriy) = 



Ittkso 



w/2 

/ 



8c7-(j'') In 



y -y 



dy' , (10) 



2 - Jmax = 0.441 (xm. 



(9) 



Thus, a rather large 0.81 (xA current circulates around 
the sample at 12.3 T when /sd = |xA, A = 0.5 (xm and 
Vm = 0.75 V. The maximum extent of this current is 60 
times farther from the sides of the sample than that 
produced by skipping orbits bouncing off of a hard wall 
with a cyclotron radius or magnet length Is = {hleUy^ 
of 7.3 nm. 



where 



haiy) 



'hB 



d/ 



V^iy) 



le 
hcOc 



d/ 



Vr(y), (11) 



as shown in Appendix A. Riess [20] extended this po- 
tential to a 2DEG with finite thickness. Thouless [21] 
then found an analytic logarithmic approximation of this 
potential far from the sample sides, and Beenakker and 
van Houten [22] then approximated the near-edge be- 
havior by introducing a cut-off at a distance ^ from the 
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sample side, and a linear extrapolation for 



Vriy) 



w/2 

/ 



Ay" 



V.iy') In 



w 



\y-y\ 



Ay (12) 



from \y\ = wl2- ^io\y\ = wl2. The characteristic length 
^is ^= UIIttu* = ie^KlTTKeJiw^ ^ where Ib = (hleB'f'^ 
is the magnetic length and a* = A7TKSaff'lm*e^ is the 
effective Bohr radius in SI units. Our values of ^, Zb , 
and a* for the ( = 2 plateau at 12.3 T are 3.3 nm, 7.3 nm, 
and 10.2 nm, respectively. 

Balaban, Meirav, and Shtrikman [10] used a nonlin- 
ear (quadratic) extrapolation near the sample sides and 
obtained the charge-redistribution potential 



^(j) 



IsdRh 



In 



w 



2^ . 



In 



y + w/2 
y - w/2 



(13) 



for ly I < w/2 - S, where 8 = Ib for the ( = 2 plateau, and 
8 is not the differential 8 of Eq. (10). They successfully 
used this potential to describe the sample-width depen- 
dence for breakdown at small currents, but could not 
account for the larger breakdown currents observed in 
other experiments [11, 23-28]. Their geometry factor is 



\ w ^+8 
In ^ -H ^^^ 



0.08 



(14) 



8 2^ 
for our values of ^ and S at w = 400 (xm. 

4.1 Charge-Redistribution Potential Equation 

The charge-redistribution potential described by 
Eq. (13) was calculated for an infinite square- well con- 
fining potential, and must be modified for use with a 
more realistic confining potential. To do this correctly 
would require a numerical solution of Eq. (12), with the 
confining potential included, as is discussed in 
Appendix A. We approximated this numerical solution 
(and then tested the approximation) by using the form of 
the potential in Eq. (13) but introducing two parameters, 
jniin and jniax, that alter the charge-redistribution poten- 
tial due to the presence of the quadratic confining 
potential. 

It was necessary to do this because the potential 
distribution of Eq. (13), with a cut-off distance 8 = Ib, 
gave the correct quantum Hall voltage Vh = Rnho across 
the sample, but the electric field E^ =- VVr did not 
increase quickly enough for increasing current to satisfy 



the /sD = (xA conditions of Sec. 3.3 and then reach 
the electric field values necessary for quasi-elastic inter- 
Landau level scattering (QUILLS) transitions [11, 16- 
18,25-28] at high currents. 

We use the same form for the charge-redistribution 
potential as Balaban et al. [10], but with a different 
geometrical factor and very different cut-off values, ymin 
and jniax, which vary with applied current. Our charge- 
redistribution potential is 



Vr(y) = -^ 



IrRn 

2 



In 



-H w/2 



w/2 



for 



where 



w 



<yn, 



■yrr 



y^yrr 



In 



y + w/2 
y -w/2 



w 



(15) 



/r = /sD - /c (rhs) - h (Ihs). 



(16) 



/c(rhs) and /c(lhs) are defined by Eqs. (5) and (6), and 
the geometry factor G in Eq. (15) is 



G(W, Jniax) = 



In 



Jmax + w/2 

w/2- y^^ 



(17) 



We assume G is current-independent, and assign the 
value 



0.147 



(18) 



to Eq. (17) by using the value of y^ax = 199.559 (xm 
found in Sec. 3.3 for /sd = (xA and w = 400 (xm. 
Our value of G is thus somewhat larger than the value 
G = 0.08 that would be used by Balaban et al. [10]. The 
cut-off values 



w/2 -y max and Sn,i„ = w/2 + y^ 



(19) 



will be determined in Sec. 5. Appendix B discusses the 
agreement between our Eq. (15) and the self-consistent 
Eqs. (10) and (11). 

The electric field E^ = -\Vr due to redistribution of 
the 2DEG with applied current is 



w 



^^^)-V^[(w/2)^-y^] 



(20) 



We now have nearly all the information necessary to 
determine the potential and current distributions. 
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5. Calculations 

Figure 4 shows the confining potential -V^(y) and 
the charge-redistribution potential -Vr(y) across the 
sample for greatly exaggerated values of A, S^ax, and 
Smin, and for an arbitrary value of /„ where /^ is defined 
by Eq. (16). Vr becomes infinite at ±w/2, but that is of 
no concern because it is only the occupied states which 
contribute to the Hall voltage, and those states occur 
only between j^ax and y^i„. The potentials are therefore 
finite and well-behaved in the region of interest. 




Fig. 4. The confining potential -Vc(y) and the charge-redistribution 
potential - Vr(y) across the sample for greatly exaggerated values of 

A S„„ = wl2 — jimnx, and S„i„ = w/2 + y^i„. 



of jmax for a particular type of experiment, and we 
believe that the results are representative of most other 
experiments since our calculations agree with the ex- 
perimental data of Fontein et al. [8]. We first note that 
E{y) = - VV{y). Therefore 



St (y max) — ^cCymax) + ^rCymax) 



(23) 



In an experiment described in Ref. [11] we measured 
the quantized longitudinal voltage drops along a GaAs/ 
AlGaAs sample between potential probes 4 and 6 of 
Fig. 1 at high currents, and deduced the maximum elec- 
tric field iimax from a quasi-elastic inter-Landau level 
scattering model. The results were 



1.1X10*^ V/m @ /sD = 215 (juA (24a) 



and 



£„ax = 4.2X10*^ V/m @ /sD = 225 (xA. 



(24b) 



The value Emax = 1.1X10'' V/m at /so = 215 [xA was just 
sufficient to excite the lowest, M = 1, QUILLS transi- 
tions [11,25-28]. It is clear from Fig. 4 that E^ax will 
occur at j^ax, so 



^tO^maxJ — ^m 



(25) 



5.1 Total Potential 

Of course the electrical transport properties depend 
on the total potential Vt(y), but we can unambiguously 
separate Vfy) into the confining and charge-redistribu- 
tion potential components 



We can therefore use Eqs. (23) and (24) to determine 
jmax, and then Eq. (22) to obtain j^in for the sample of 
Ref. [11]. Note that changing the values of j^ax and j^in 
also alters the values of /^(rhs), /c(lhs), and thereby the 
value of /r in Eqs. (5), (6), and (16). Thus there are no 
free parameters , and one can obtain unique solutions to 
the total potential and other transport properties. 



Vt(y) = Vc(y) + Vr(y) 



(21) 



5.2 Results 



We have uniquely defined the potentials V^iy) and VXy) 
in Eqs. (2) and (3) of Sec. 3.2 and Eqs. (15) to (18) in 
Sec. 4.1, plus Eqs. (5) and (6) in Sec. 3.3. The current- 
independent parameters for the confining p otential and 
the charge-redistribution potential are: A = 0.5 [xm, 
Vm = 0.75 V, and G = 0. 1 47. For a given sample we know 
the applied current /so and the sample width w , but there 
are still two free parameters: j^ax and y^^- 

Ordinarily, it would not be possible to uniquely deter- 
mine the values of j^ax and y^^ since the only other piece 
of information is that the quantum Hall voltage Vh is 



Vh = ^h/sD = Vt(ymi„) - VtCymax), 



(22) 



and there is a range of values for jmax that satisfies this 
equation. It is possible, however, to determine the value 



Relevant values for the solution at /so = [xA are 
shown in Table 1. Most were calculated in Sec. 3.3; the 
remainder were found from Eqs. (2), (3), and (15) to 
(22). Note that j^ax and y^^ are predicted to be about 60 
magnetic lengths from the sides of the sample. 

We calculate the values shown in Table 1 at 
/sD = 215 [juA by increasing the value of j^ax until 
£t(ymax) = 1.1 X 10'' V/m, adjusting the value of jmin to 
obtain the correct Hall voltage, and remembering that 
changing the values of j^ax and j^in also changes the 
values of /^(rhs), /c(lhs), and I^. The solution is unique, 
with no free parameters. The same procedure is done at 
/sD = 225 [juA, except that the value of j^ax is increased 
until EtCymax) = 4.2 X 10'^ V/m. Note in Table 1 that y^„ 
is still about 13 magnetic lengths away from the side of 
the sample at /so = 225 [xA. 
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Table 1. Values of some quantities obtained from the calculations in Sec. 5 for Isd = [lA, 25 (xA, 215 (xA, and 225 (jlA. The quantities 
common to all the calculations are / = 2, B = 12.3 T, k = 13.1, w = 400 (xm, A = 0.5 (xm, y„ = 0.75 V, a = 3.0X lO'^ VW, A = 199.500 
(xm, G = 0.147, and Is = 7.3 nm. See Sees. 2-5 for the definitions of these symbols 





/,(rhs) 


/,(lhs) 
(t^A) 


(H,A) 




Jmin 


(V) 


(V) 


£c(y™x) 

(MV/m) 


£r(y„ax) 

(MV/m) 


£.(y™x) 

(MV/m) 


S^JIb 





0.81 


-0.81 


0.00 


199.559 


-199.559 


-0.010 


0.000 


0.354 


0.000 


0.354 


60.3 


25 


0.94 


-0.68 


24.74 


199.564 


-199.554 


-0.012 


-0.160 


0.382 


0.054 


0.436 


59.6 


215 


2.30 


-0.05 


212.75 


199.599 


-199.515 


-0.030 


-1.392 


0.596 


0.504 


1.100 


54.8 


225 


37.36 


-0.00 


187.64 


199.901 


-198.044 


-0.482 


-1.477 


2.405 


1.795 


4.200 


13.5 



We also calculate the relevant quantities at /so = 25 
(jlA, which is a current often used in precision quantized 
Hall resistance measurements. In this case, however, we 
do not know the value of iitCymax), so we use a linear 
interpolation of the value of y^ax between its values for 
/sD = (juA and 215 (xA. The quantities shown in Table 1 
for /sD = 25 (juA are relatively insensitive to this choice 

for Jniax- 

5.3 Plots 

We now plot the potentials, using Eqs. (2), (3), (15) 
to (19), and (21). Figure 5 shows V^(y) and Vr(y) for the 
parameters used in Table 1 at /so = 215 (juA, except that 
the plot is between ±0.99999 w/2(± 199.998 (jum) 
rather than y^ax and ymin in order to show the sharpness 
of the confining potential and the extent of the charge- 
redistribution potential at these extreme values of y. 
Figure 6 shows Vt(y) plotted between y^ax and ymin using 
the parameters in Table 1 at /jd = 215 (juA and 225 (juA. 
Other than moving farther to the right, the total potential 
does not significantly change shape with increasing 
current. 

Figure 7 shows Vt(y)at /so = 25 (juA. The shape of this 
predicted potential is in excellent agreement with the 
experimental measurements shown in Fig. 6 of Fontein 
et al. [8]. It is this agreement which provides the best 
verification of our results. The "linear" part of the po- 
tential distribution within the sample interior, attributed 
in Ref. [8] to heating effects which cause R^ = VJho to 
increase, is accounted for by our charge-redistribution 
potential in a sample which has minimal heating at these 
currents [24]. 

The electric fields E^iy) = - VVcCy) and E^iy) = - 
VVrCy) are shown in Fig. 8 for /jd = 215 (juA; they were 
determined from Eqs. (2), (3), (18), and (20). The value 
of y:nax= 199.599 |xm is such that Efy^,^) = \ .\X\(f 
V/m in equation (23). The contribution to the total elec- 
tric field at y^ax is slightly more for the confining 
potential than for the charge-redistribution potential at 



this current. Table 1 shows that the confining potential 
also provides the dominant contribution to StCymax) at 
other currents. 




V,(y) @ 215 nA 



V^(y) @ 215 fiA 



-200 



-100 



100 



200 



y (lim) 



Fig. 5. K(y) and V,(y) plotted between ±0.99999 w/2 for the 
parameters used in Table 1 at ho = 215 (jlA. The parameters common 
to all plots in Figs. 5-1 1 are / = 2 (12 906.4 fl), B = 12.3 T, w = 400 
(jLm, K = 13.1, A = 0.5 n,m, y„ = 0.75 V, and G = 0.147. 



The location, y^ax, of the last- filled state on the right- 
hand side of the sample increases with applied current 
/sD. We can use Eq. (A-3) and Table 1 to determine what 
part of this increase in y^ax is due to the increase in the 
total electric field at y^ax- The percentage contributions, 
relative to the values of y^ax and iitCymax) at /so = (juA, 
are 4 %, 5 %, and 3 %, for ho = 25 [lA, 215 [lA, and 
225 (xA, respectively. Therefore, most of the increase in 
y^ax is due to the Lorentz force pushing the electrons 
closer to the side of the sample. 

The current density Jt(y) for electrons moving in the 
positive X direction is 



/.(y) = a^Efy) = — [£,(y) + E,(y)] 



(26) 
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vitxv© aasM 




Fig. 6. ViCv) plotted between y„ax and y„i„, using the parameters in 
Table 1 for /sd = 215 (jlA and 225 \xA. The values of y„„x and y„i„ are 
199.599 (JLm and -199.5 15 (xm, and 199.901 (jlui and -198.044 (xm for 
/sD = 215 (xA and 225 [lA, respectively. 
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Fig. 7. Vi(y) at /sd = 25 (jlA. This potential is in excellent agreement 
with the experimental measurements shown in Fig. 6 of Fontein et al. 
[8]. 



Figure 9 shows Jt(y) for Igo = 25 [xA, 215 [xA, and 225 
[xA. The maximum two-dimensional current density is 
at jniax, and is 85 A/m and 325 A/m, respectively at /sd 
= 215 [xA and 225 [xA. There is current in the negative 
X direction in the vicinity of jmm at small currents due to 
the dominance of the confining potential. When 
/sd = 215 (xA and 225 (xA, however, /irCymm) > \E^(y„,a)\ 
and no current flows in the -x direction anywhere across 
the sample. 
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Fig. 9. Current density Ji(y) for hu = 25 (xA, 215 (jlA, and 225 (xA. 
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Fig. 8. Electric fields fi^Cv) = - VK(y) and £r(y) = - Vy,(y) for 
/sd = 215 n,A. 



The current /(jy) for electrons moving in the positive 
X direction is 



}/.( 



ny) 



i(y) = j J^(y)dy = -^, (27) 





where 



/sd = J /tO')d>' = /(J^ax) + liymin) , (28) 

Jmin 



and 



M = I(y2)-I(yi). 



(29) 
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We divide the sample width into 20 equal segments in 
Fig. 10 and determine the percentage of current flowing 
through each segment for Isd = 25 [xA, 215 [xA, and 225 
[xA. We do not show a plot for Isd = [xA, but A/ would 
be -0.81 fxA and +0.81 fxA for the left-hand side and 
right-hand side segments, respectively, and zero for the 
other 18 segments because /c(rhs) = -/c(lhs) = 0.81 (xA. 



Finally, we investigate the charge-redistribution 
-e^aiy) of the electrons in the 2DEG in terms of the 
deviation 8(T(y) in the number of electrons/cm^ from 
the average number ris = 5.94X10'Vcm^ on the ; = 2 
plateau at 12.3 T, where 






(30) 
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Fig. 10. Percentage of current flowing through 20 equal segments 
across the sample width for Isd = 25 (xA, 215 (xA, and 225 [lA. 



from Eq. (A-5). Figure 11 is a logarithmic plot of 
I8cr(3;)| versus y for /so = 215 [xA. There is an excess of 
electrons on the +y side of the sample, and a depletion 
on the -y side. 

An assumption made in deriving Eq. (A-5) was that 
the charge density varies slowly across the sample, i.e., 
that p(y) ~ p(y+dy), or 8(T(y) « 11^. This assumption 
is valid here because the largest value of ^cr{y) occurs 
at jmax , and is 2 %, 2 %, and 6% of n^ at ho = 25 [xA, 
215 (xA, and 225 [xA, respectively. 

One of the consequences of our approximate form of 
the charge-redistribution potential is that the net charge 
does not vanish when the charge-redistribution -e^cr{y) 
is integrated across the sample width. The area under 
the curves in Fig. 11 is 4 % larger for the +y side than 
for the -y side. Therefore, there is an unaccounted 
excess of electrons; so this is not quite the actual shape 
of the charge-redistribution function. However, it is the 
potential and current distributions that are of primary 
importance to the transport properties — not the charge- 
redistribution. The charge was certainly conserved in 
the experiment of Fontein et al. [8], and yet their mea- 
sured potential distributions are symmetrical. This fact 
demonstrates that the slight charge asymmetry does not 
significantly affect the potential and current distribu- 
tions. 



The current distributions in Fig. 10 are virtually iden- 
tical between 25 [xA and 215 (xA, even though large 
numbers of electrons are being excited into higher 
Landau levels at 215 [xA. The left and right side distri- 
butions are nearly symmetric. There is, however, a sig- 
nificant transfer of current from the left-hand side seg- 
ment to the right-hand side segment at 225 (xA. We saw 
in Sec. 5.2 that no current flows within 60, 55, and 13 
magnetic lengths of the sample side for ho = 25 [xA, 
215 fxA, and 225 |xA, respectively Also, 68 %, 70 %, 
and 51 % of the current is in the 19 segments to the left 
of the right-hand side segment where the edge channel 
current would flow for these three applied currents. The 
current density was negative in the left-hand side of 
Fig. 9 at Isd = 25 [xA because electrons were flowing in 
the -X direction at y^^ -199.554 [xm, but that contribu- 
tion to A/ in the left-hand side segment of Fig. 10 is so 
small that the net current is positive. 




Fig. 11. Logarithmic plot of the charge-redistribution IBtrC^;)! across 
the sample for ho = 215 (xA, where 8(7(j') is the deviation from the 
average number density n^. The "-" region represents an excess of 
electrons, the "+" region a depletion of electrons. 
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We could conserve the charge by adjusting the origin 
slightly to the right until the area under the curves are 
equal for ±ym Fig. 11, and then self-consistently recal- 
culating the potentials with the new coordinates. This 
would greatly complicate the calculations however, and 
with all the approximations that have been made in this 
paper, and with the excellent agreement with experiment 
[8], it seems unnecessary. It may be a consequence of 
this charge nonconservation problem that the value of 
jniin is inside the confining potential for the case in 
Table 1 when /sd = 225 (x A. 



5.4 Sample-Width Dependence 
Current 



of the Critical 



Balaban, Meirav, and Shtrikman [10] have found that 
the critical current for breakdown of the quantum Hall 
effect, /cr, scales logarithmically with the sample width 
w for all Landau levels. We verify this dependence by: 

(a) using the result in Sec. 3.3 that /c(rhs) = - /c(lhs) = 
0.81 (xA for the i = 2 plateau at 12.3 T when /sd = (xA; 

(b) calculating the value of j^ax from Eqs. (2) and (5) for 
each value of w at /sd = (xA; (c) calculating the value 
of G from Eq. (17) for each value of w ; (d) defining /cr 
as the applied current /sd sufficient to excite the lowest, 
M = \, QUILLS transitions [11,25-28]; (e) assuming 
the value of /ScCVmax) is the same for all values of /sd that 
excite M = \ QUILLS transitions (The value used is 
ficCVmax) = 5.96X10^ V/m, obtained from Table 1 at 
/sd = 21 5 (xA, w = 400 (xm, and E^iy^^) =1.1X10" V/m. 
This is equivalent to fixing the value of /^(rhs) to be 2.30 
(xA for each value ofI„); (f) calculating the value of j^ax 
from Eq. (2) for each value of w ; (g) adjusting the value 
of/, so that EtCVmax) =1.1X10" V/m in Eqs. (2), (17), 
(20), and (23); and (h) adjusting the value ofy^m to give 
the correct Hall voltage for each current by using Eqs. 
(2), (15), (17), (21), and (22). 

The results of I„ versus w are plotted in Fig. 12. The 
shape of the curve is identical to the experimental data 
of Balaban et al. [10]. The scaling is very different, 
however because their critical currents are about two 
orders of magnitude smaller than ours. We note that the 
experiment of Haug, von Klitzing, and Plog [29] tends 
to agree with the experimental curve shapes of Balaban 
et al. [10], but the experiment of Kawaji, Hirakawa, and 
Nagata [30] found a linear, rather than a logoarithmic, 
dependence of I„ with w . Perhaps this difference is due 
to nonuniformities in the values of the charge-depletion 
width A along the sides of the samples, e.g., we have 
observed different values ofI„ along the lengths of some 
of our samples. If we assume that the value of V^ 
remains constant along a sample edge, allow A to vary 
by changing the average ionized donor density A^d, and 
assume the ratio £'cCyraax)//itCyraax) remains constant, then 



we find that the critical current required to excite M =\ 
QUILLS transitions with £,Cyn,ax) = 1.1 X 10" V/m 
decreases when A decreases and A^d increases, i.e., the 
steeper the confining potential, the smaller the critical 
current. 
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Fig. 12. The critical current /„ versus sample width w. Refer to Sec. 
5.4 for details of the calculation. 



6. Conclusions 

We have calculated potential and current distributions 
across the width of a GaAs/AlGaAs heterostructure 
sample for applied currents between (xA and 225 (xA, 
using: (a) a quadratic confining potential KCv) arising 
from charge-depletion regions along the sides of the 
sample; (b) parameters for that potential obtained from 
a localization experiment [9]; (c) a logarithmic charge- 
redistribution potential V^iy) of the 2DEG; and (d) a 
maximum electric field /StCVmax) calculated from break- 
down measurements and a QUILLS model [11]. Our 
predictions are in excellent agreement with experiments 
[8,10]. 

Referring to Table 1 , the confining potential compo- 
nent /icCVraax) of the electric field at j^ax contributes 
88 %, 54 %, and 57 % to EtCVmax) at 25 |xA, 215 (xA, and 
225 (xA, respectively. The maximum current density 
•/tCVraax) is 34 A/m, 85 A/m, and 325 A/m, respectively at 
these three currents. A significant amount of current is 
distributed within the sample interior. For example, /, is 
99 %, 99 %, and 83 % of /sd, respectively at these three 
currents. We predict the current to be much farther from 
the sides of the sample than in other models, e.g., no 
current flows within 60, 55, and 13 magnetic lengths of 
the sample side for these currents. It would require a 
lateral resolution of about 0. 1 (xm, rather than the 70 (xm 
resolution of Fontein et al. [8], to verify this result. 
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7. Appendix A. Derivation of Eqs. (10) 
and (11) 

We first derive Eq. (10) of Sec. 4. The Lorentz force 
exerted on the conducting electrons causes deviations 
-eh<j(y') from the average surface charge density 
-eris = - ie ^B/h of the 2DEG at each point y ' across the 
sample width. Consider a strip of this redistributed 
charge of width dy', located in the x-y plane at position 
y' and pointing in the x direction, with a charge/length 
A (y') = -eh<j(y')dy'. Gauss's law is then used to obtain 
the electric field dE(y) at some point y in the x-y plane, 



due to one of these line charges: KSo'fE • d^ = q, and 

thus KSodE(y)2TT\y'-y\Lx = A(y')Lj, for a cylindrical 
Gaussian surface of radius \y'-y\ and length L^. The 
potential dV(y) of this line charge is 



dV(y) 



-/' 



dE(y)dy = 



A(y') 
2tj-kso 



w J 



1 



y'—w/2 



2\y'-y\/w 



dy 
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Fig. A-1. Geometry factors for the calculations in Appendix A. 

If the volume charge density is slowly varying, then 
p(y) ~ p(y + dy), where ph = -ecr and 8p/i = - eScr. 
Thus, from Eqs. (A-3) and (A-4), 

t>CT{y)dy = ^ [E{y)-E{y + dy)], 



. ,,_ ajy) dE(y) ajy) d'V(y) 
^'^^^- oj,B dy ~ w,B dy^ 



The total potential, when summed over all the line 
charges, is Eq. (10) 



V{y) = 



Now we derive Eq. (1 1) of Sec. 4. We found in Eq. (3) 
of Ref . [ 1 1 ] that the center of mass coordinate yo of each 
state undergoing cycloidal motion in a Landau level is 



e 


wl2 

8cr(y') In 

-w/2 


- \y'-y\ 

_w _ 


dy. 

(A-2) 


2t7KSo 



_ E(y) M^ _ 
^° wM eB ' 



(A-3) 



so the states move to the right as ^(y) increases. The 
total charge SQ transferred into the volume l^hdy, out- 
lined with solid lines in Fig. A-1, is 



hQ = hphhdy = lMAyo(y)p(y) 



-Ayo(y + dy)p(y + dy)]. 



(A-4) 



But <T(y) = ieB/h, so we obtain Eq. (11) 
ie d^Vjy) _ im* d^Vjy) 



8(7 (y) 



hu>^ dy 



hB dy^ 



(A-5) 



Equation (11) considers the charge-redistribution 
-e8cr(y) due to the second deriwtive of the charge- 
redistribution potential Vriy). The charge-redistribution 
that we calculated in Sec. 5.3 depends also on the 
second derivative of those regions of the confining 
potential V^{y) which differ from the y^ax and y^,^ values 
at ho = (xA. 

8. Appendix B. Eigenvalue Equation 

We saw in Eq. (12) of Sec. 4 that 

"2 



Vr(y) = -^ j ±r2 [Vriy')] In 



dy' 



w 



\y'-y\ 



where the characteristic length ^ is 



2TTKSohcO^ 



dy', (B-1) 



(B-2) 
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MacDonald, Rice, and Brinkman [19] have pointed out 
that Eq. (B-1) has the form of an eigenvalue equation 

@[F(y)] = CFiy), where 



W = - J^[]ln 



\y'-y\ 



Ay' 



is an operator acting on the function F(y) = Vr (y'), and 
C = 1/^ is a constant. 

The function F(y) = Vr(y') should have a second 
derivative that satisfies Eq. (B-1). Our function, from 
Eq. (15) of Sec. 4.1, is 



Vr(y') = -^-f^ 



In- 



In 



y' + w/2 
y' - w/2 



, (B-3) 



where the cut-off is S^ 



w 



d 17 ^ IN ItRh 



In 



w 



[iw/2r-(yr 



and 



d T7 / i\ ItRh 



In 



2wy' 



[iw/2f-(yrf 
(B-4) 



We can see if the potential given by Eq. (B-3) is a 
valid solution to the eigenvalue expression Eq. (B-1) by 
substituting Eq. (B-4) into Eq. (B-1) and integrating 
only between the limits - S^in to S^ax because 8(t()'') is 
zero beyond these two cut-off values. Surprisingly, we 
obtain nearly exact solutions to the eigenvalue ex- 
pression at all values ofy when ^is less than 2.0 |jum and 

Smax = Smin = ^- ThlS CholCC of Sn,ax = S^in = ^ fOT thC 

cut-off values was used by Beenakker and van Houten 
[22]. Our value of ^ is 3.3 nm, so we would be well 
within this exact range if Sn,ax = S^in = ^. The values of 
Smax do not equal Sn,;n however, and are much larger than 
the value of ^. Also, we have a confining potential, 
Vc(y), parts of which should be included in Eq. (B-1). 
It would be interesting to see how well the values of 
Vt(y) obtained in Eq. (B-4) agree with the values 
obtained in the eigenvalue expression Eq. (B-1) when 
using the values of quantities obtained in Table 1. 
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